A theory of the three-pulse electric-dipole echo in glasses in a magnetic field 
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We extended existing theory of the two-pulse electric-dipole echo in glasses in a magnetic field 
to the three-pulse echo. As is well known, at low temperatures two-level systems (TLS's) are 
responsible for the echo phenomenon in glasses. Using a diagram technique in the framework of 
perturbation theory we derived a simple formula for the three-pulse echo amplitude. As in the case of 
two-pulse echo the magnetic field dependence of the tree-pulse echo amplitude in glasses is related 
to quadrupole electric moments of TLS's non-spherical nuclei and/or dipole-dipole interaction of 
their nuclear spins. These two mechanisms are responsible for the fine level splitting of TLS. As a 
result TLS transforms to multi-level system with the fine level splitting depending on a magnetic 
field. Due to existence in the theory the additional parameter T — the time interval between the 
second and the third pulses we have more reach spectrum of echo oscillations in a magnetic field in 
comparison with the case of the two-pulse echo. 

PACS numbers: 61.43.Fs, 76.60.Gv, 77.22.Ch 
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I. INTRODUCTION 

It is well known that glasses at low temperatures (be- 
low 1 K) display a number of universal properties which 
are drastically different from the properties of similar 
crystals. These properties are almost independent of the 
chemical composition of a glass and mainly due to the 
disordered atomic structure of an amorphous solid; more 
exactly, due to the absence of a long-range order in the 
glass 

According to the standard theory of Anderson, 
Halperin, Varma, and Phillips all these universal 
properties are associated with the existence of two-level 
systems (TLS's) in glasses. TLS consist of atoms or 
groups of atoms that may tunnel between two minima in 
a double- well potential. However, the microscopic nature 
of the TLS's in the majority of glasses remains unknown. 

One of such universal phenomena is an electric-dipole 
echo [l|. Depending on the sequence of applied high fre- 
quency electric pulses one usually discriminates between 
two-pulse and three-pulse echo. The essence of the two- 
pulse echo is the following. When a glass is subjected to 
two short electromagnetic pulses at a frequency of about 
1 GHz separated by a time interval t, a response signal 
at the same frequency may be detected in the dipole mo- 
ment of the system at time r after the second pulse. In 
the three-pulse echo we have a sequence of three pumping 
electromagnetic pulses in the moments 0, r and t+T cor- 
respondingly. A coherent echo signal in this case appears 
at the time T + 2r. 

A very interesting echo phenomenon was discovered re- 
cently in 2002 [2j . The amplitude of the two-pulse dipole 
echo in some nonmagnetic glasses exhibited a strong non- 
monotonic (oscillating) dependence on the magnetic field 
even in weak fields (about 10 mT). Wurger, Fleischmann, 
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and Enss Q suggested that this unusual magnetic field 
effect is caused by the presence of non-spherical tunneling 
nuclei with the electric quadrupole moment (or, equiva- 
lently, the nuclear spin J > 1) in the glass. The inter- 
action of the quadrupole moment with gradient of static 
microscopic electric field in two wells of double-well po- 
tential creates a fine splitting of energy levels in TLS 
transforming the TLS to multi-level system. 

This idea agrees very well with the experimental data, 
in particular, with the results of the recent measure- 
ments of the two-pulse dipole echo in glycerol (CsHsOa) 
4]. In these experiments, it was shown that the replace- 
ment of hydrogen, which has nuclear spin J — 1/2 and, 
consequently, a zero quadrupole moment, by deuterium 
(J = 1 and a nonzero quadrupole moment) results in 
a more than an order of magnitude enhancement of the 
magnetic-field dependence of the echo. 

The magnetic-field dependence of the two-pulse echo 
amplitude in glasses with nuclear quadrupole moments 
was theoretically studied by Wurger, Fleischmann, and 
Enss [l, [11 and by one of us [1, [3|- However the weak 
magnetic field non monotonic dependence of echo in non- 
deuterated glycerol could not be explained in such a way. 
Neither of nuclei forming glycerol glass CsHgOa possess 
a nuclear quadrupole moment. 

To explain this interesting phenomenon in non- 
deuterated glycerol Bazrafshan et al. [Sj suggested that 
the magnetic field dependence of the two-pulse echo am- 
plitude in this case may be caused by the dipole-dipole 
interaction of the nuclear magnetic moments of hydrogen 
atoms. This interaction (along with Zeeman interaction) 
can also create a fine structure of two levels in TLS, which 
depends on the applied magnetic field. 

In the papers [1, the two-pulse echo amplitude was 
numerically calculated in deuterated glycerol C3D5H3O3 
(taking into account all of the spins of hydrogen and the 
quadrupole moments of deuterium nuclei), assuming that 
the tunneling motion in the two level system is the ro- 
tation of the glycerol molecule as a whole. Results of 
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FIG. 1: Electric pulses and the three-pulse echo in the TLS. 



these calculations were in a good agreement with exper- 
imental data. Later in an analytical theory of the 
magnetic field dependence of the two-pulse echo ampli- 
tude in glasses with dipole-dipole interaction of nuclear 
spins was developed. Without any assumptions about 
microscopic tunneling mechanism of hydrogen atoms in 
glycerol and without fitting parameters it also shows a 
very good agreement with experimental data. 

The purpose of this paper is to build a theory of the 
three-pulse echo in glasses in a magnetic field to stimu- 
late further experimental investigations of this interesting 
phenomenon. We consider here two mechanisms of mag- 
netic field effects, namely quadrupole electric moments 
of TLS's nuclei and dipole-dipole magnetic interaction of 
their nuclear spins. 



dipole moment of the TLS, Aq is the tunneling splitting, 
E = E2 — El is the TLS energy, and Ei and E2 are the 
energies of the ground and the excited states of the TLS, 
respectively. The electric field frequency is assumed as 
huj « E. To simplify our equations we will put h — \. 

After the action of the 1st electric pulse C2 acquires 
in the first approximation a finite value proportional to 
amplitude Vi which is assumed to be small. During the 
time interval between the first and the second pulses (0 < 
t < r) we have 



Ci oc 1 • e 



C2 OC -iVie-'^^\ 



(5) 



where we have neglected terms proportional to V-^. In a 
similar way after the action of the second pulse we have 
between the second and the third pulses 



C2 cx -iVie-^^^^ - iV2e 



(6) 



Then finally after the third pulse we get 



-iE^iT+T) -iEi{t-T-T) _ 
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II. THREE-PULSE ECHO IN AN ENSEMBLE 
OF TLS 

We start with a simple example of a three-pulse echo in 
an ensemble of TLS. There are three exciting electromag- 
netic pulses with amplitudes Fi, F2, and F3 correspond- 
ingly. The time interval between the first and the second 
pulse is r, between the second and the third pulse — is 
T (Fig. [1]) . We suppose that the duration of any pulse is 
much less than the time delays between the pulses. 

First we consider one TLS. The wave function of the 
two level-system is a linear combination of wave functions 
Lpi and (^2 for each level 



^ ^ Ciipi + C2ip2, \Ci\ 



\C2? = l. 



Initial conditions (just before the first pulse) are 
Ci = 1, C2 = 0. 



(1) 



(2) 



Then in the electric field the time variation of the ampli- 
tudes Ci and C2 obeys the equations 

ih'^ = EiCi + v{t)C2, ih^ = E2C2 + v{t)Ci. 

at at 

(3) 

For the off-diagonal transition matrix element V{t) we 
have a following expression during the electric pulse 

V{t) = Fi,2,3 costot, Fi,2,3 = (Fi,2,3 ' m) ^ . (4) 

Here Fi^2.3 is the electric field amplitude of the first, sec- 
ond and the third electric pulse, respectively, m is the 



C2 cx ^tVie-'^^' - zl/2e-*^^"e-'^^(*-")- 

-iy3e-*^i(^+^)e-*^^(*"^-^). (8) 

The complex echo amplitude is proportional to C1C2. 
Only a product of the second term in Eq. ^ and the 
third term in Eq. ^ and a product of the third term 
in Eq. ([7]) and the second term in Eq. ^ gives us the 
three-pulse echo signal. As a result of these two contri- 
butions the three-pulse echo amplitude from one TLS is 
proportional to 



(9) 



Summing echo amplitudes over all TLS's we have for the 
total echo amplitude 

Pocho cx ^ ViV2V,e^^^'-^-'^\ (10) 

TLS 

At the time t = T + 2t the phase of the echo signal from 
all TLS's becomes equal to zero so we have a maximum 
of the overall echo amplitude Pecho cx —iViV2V3N, where 
N is the number of resonant TLS's (i.e. TLS's with hu « 
E). 



HI. DIAGRAM REPRESENTATION OF THE 
ECHO 

For the future analysis of the echo in magnetic field let 
us now give a simple diagram representation of the three- 
pulse echo from one TLS (see Fig. [J). For simplicity, we 
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FIG. 2: Diagrams representation of the three-pulse echo from a TLS. (a) The first contribution, (b) The second contribution. 



assume that the electric field amplitudes Fi^2,3 of the 
three electric pulses are parallel to each other, i.e. = 
Fie, where e is a unit polarization vector. 

Diagram in Fig. [2K corresponds to the product of the 
second term in Eq. ^ (the full line in the diagram) and 
the third term in Eq. ([5]) (the dashed line). Each vertical 
line (except the final wavy line) corresponds to a transi- 
tion from one level to another under the action of one of 
electric pulses. One should ascribe to this line a factor 
{—iFj)aku where j is a number of the pulse, and au is 
a dipole transition matrix element from the level / to the 
level k in the TLS. Each horizontal line corresponds to 
free TLS evolution in the time intervals between or after 
the pulses. For example between the neighboring pulses 
i and j we ascribe to such line a factor ex.p{—iEkStij), 
where A: is a level number, and dtij is a time delay be- 
tween the i and j pulses. Finally factor a2i corresponds 
to the vertical wavy line. In the final expression all fac- 
tors corresponding to dashed lines should be taken as 
complex conjugated [S]. 

According to these rules we get for the first contribu- 
tion to the echo amplitude the following expression 



{-iF2)ai2 ■ e 



a2i 



(11) 



We took here into account that ai2 = a*2x The second 
contribution corresponds to the product of the third term 
in Eq. ([7]) and the second term in Eq. ([5]) 



xe 



'iEi{t-T-T) 



^^(^+^).(-*F3)ai2X 



^21 



^iFiF2F3|Q;i2| |a2i| e 



2 i(E2-Ei)(t^T^2T) 



(12) 



Comparing Eqs. (fTTI) and (fT^ we see that both con- 
tributions to the echo amplitude turn out to be identical 
in this simple case. As a result the total echo amplitude 



Pccho oc -iFiF2F3|ai2| |a2i| e 



2 i{E2-Ei){t~T-2T) 



(13) 



IV. ECHO IN A MULTI-LEVEL SYSTEM 

Using these simple diagram rules we can now easily 
find the contributions to the three-pulse echo signal from 



a multi-level system. Let us, for example, consider a 
multi-level system as shown in Fig. |31 It consists of two 
identical groups of Af levels shifted vertically against each 
other by some energy E (playing the role of the usual TLS 
energy). Inside the groups the positions of the levels are 
arbitrary, i.e. they are not necessarily equidistant. As 
we will see in the next sections, such a multi-level system 
describes a TLS with a quadrupole nuclear moments or 
with dipole-dipole interaction of its nuclear spins. 

According to the diagram rules formulated in the pre- 
vious section we get two contributions to the echo ampli- 
tude 
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-iFi)a,2.ii ■ e' 



iEj2T ^ 



X (-ii^2)«fcij2-e-'^"(*-")x 
^e-^-(^+^).(-*F3)a,2,a-e-^-(*- 



-T) 



Oil2,kl 



FiF2Fj,aj2,iiaki,j2a*i2,iiOti2.ki'^ 



X exp[-i(£;j2 - Ea)T ~ i{Eki - Eii)T + 

+ i{Ei2~Eki){t~T-T)], (14) 



and 
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-iFi)a.2.ii ■ e 



'il,fclj"2,i2V^/ 

x(-*F2)a,i,,2-e-^"(*---^): 



-iEj2{T+T), 



<^l2,kl 



= --j^FiF2F^aj2,iiaki,j2a*i2,iiOti2.ki'^ 

X eyip[~i{Ej2 - Ea)T + i{Ei2 - Ej2)T + 

+ i{Ei2-Eki){t-T-T)]. (15) 



The factors l/vA/" in these formulas correspond to the 
case when the low-energy group levels are equally popu- 
lated (with a probability 1 /N) and the high-energy group 
levels are empty. This is the case for low enough temper- 
atures, kT ^ E. On the other hand, to satisfy the previ- 
ous conditions the temperature must be much larger than 
the width of the energy splitting in the groups (this width 
is of the order of the level splitting AE'), i.e. kT ^ AE. 
The latter inequality corresponds to the usual experi- 
mental situation. The two conditions are compatible if 
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FIG. 3: Diagrams representation of the three-pulse echo in multi-level system, (a) The first contribution, (b) The second 
contribution. 



E ^ AE. In the usual experiments this inequahty is 
met since the resonance frequency huj « E' 3> AE. The 
hmitation E 3> kT is not a crucial. The final result can 
be easily generalized to the case kT ~ £^ by including 
thermal occupation numbers. 

Two groups of levels are identical and only shifted by 
the energy E so we have 



Ei2 — E + Ell, Ej2 — E + Eji. 
Then we can rewrite Eqs. and as 



(16) 



(1) 



X eyi-p[iE{t - 2t - T) - i{Eji - Eii)T - 
- i{Eki - E,i)T + i{En - Eki){t - r - T)], (17) 



and 



,(2) 



lit) 



"ii,fcij2,i2VW — -j^FiF2F^a^^ ctkj^ii ^ik ' 

X eyi^[iE{t - 2t ~ T) - i{Eji ~ Eii)T + 
+ i{En - E,i)T + i{En - Eki){t - r - T)], (18) 

where we write the indexes (1, 2) of the level groups as 
superscripts. From Eqs. p?)) . p5|) it follows, in the case 
when E ^ AE, that the echo signal appears a,t t — 2t + 
T. At this time, summing over all possible combinations 
of vertical transitions (i,j,k,l) between different levels 

and taking into account that Q^ij^'' — finally get 

for two contributions to the echo signal from one multi- 
level system the following expressions 
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plcL(r + 2r)(x--.FiF2F3X 
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and 



pfcL(r + 2r)cx--^FiF2F3X 



i.k 
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We can see that these two contributions are not identical, 
as opposed to the case of simple TLS. 

Now we consider an arbitrary TLS- nuclear spins multi- 
level system. Like in [6] we use perturbation theory and 
get finally for the echo amplitude 



Pccho(T + 2t) cx -~iFiF2F3 ( ) X 



64 ^ /A 

U ^ [e 

n, ray 71 ^ 



(K 



J )nm 



sm\enmT/2h) sin2(e„„(T + T)/2h) 



(21) 



Where = Em — En is the distance between the energy 
levels inside one level group. And {Vj)nm are transition 
matrix elements which we define below. This is a general 
formula for three-pulse echo for any mechanism of TLS- 
nuclcar spin interaction. At T = and F2 ~ -F3 this 
formula transforms to the known formula for two-pulse 
echo amplitude 0]. In this sense we can say that the 
three-pulse echo is a generalization of the two-pulse echo. 



V. QUADRUPOLE INTERACTION 

In this section we specify the general formula, Eq. (|2ip , 
for a TLS with a nuclear quadrupole electric moment in 
an external magnetic field. 

The interaction between the nuclear electric 
quadrupole moment and gradient of an internal 
electric field is described by the Hamiltonian [ll[ 



E 



(1,2) 



n, (1.2) 



(1,2) _ 



^2^(1,2) 

dradrp ' 



(22) 



where 1^9^^^^) is an electrostatic potential at the nuclear 
site depending on TLS well (1) or (2). The operator of 
the nuclear quadrupole electric moment is the traceless 
tensor \\M 



(20) Qcp 



eQ 



6J(2J- 1) 



2 ( J<^Jp 



+ Jf)Jc] -5af)J{J +1) , 

(23) 
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FIG. 4: Echo amplitude as a function of the external magnetic field ^HT/2Jh for different ratios p — T/t. 



where J is the operator of the nuclear spin. Then we 
define operators 



Wj = Qcp^i+J - M • H, Vj = Qo;3^S , 

(1) _i_ (2) 
(±) _ fcp ± 
"focp - 2 ' 



(24) 



where H is the external magnetic field, M is a nuclear 
magnetic dipole moment of the TLS. 

Using unitary transformation Sj we can diagonalize 
matrix Wj and introduce a new interaction operator V j 



Sj Wj SJ^ = Wj, and Sj Vj SJ^ = Vj. 



(25) 



Now Wj is a diagonal matrix in the nuclear spin space 
Af X Af (for one nuclear spin Af — 2J + 1). ft gives us 
the nuclear quadrupole energies En and fine interlevel 
spacing = Em — En in the average (over the two 
minima) internal electric field gradient and in the exter- 
nal magnetic field H. The unitary transformation Sj can 
be found in general only by numerical diagonalization of 

the matrix Wj. The transformed matrix Vj entering 
the Eq. (PT|) describes the interaction of the TLS with 
its nuclear quadrupole moment. It is proportional to the 
difference in electric field gradients in two wells. 

Let us further consider the case when Zeeman energy 
is much larger than the quadrupole level splitting 



(26) 



and distances between the pulses satisfy the condition 



WoSiT + r)/h<l, 



(27) 



In this case we neglect the quadrupole interaction term 
Qap'f'af} (EH) in comparison with Zeeman energy. 

In a full analogy with two-pulse echo Q we get from 
Eq. (|2ip the final formula for the echo amplitude in the 



high magnetic field: 



Ao 



X U-C- 



sin^(/ii7(T -I- r)/2 J) sit? {^iHt /2J) 



sin^(/ig(r + t)/J) sm^ifiHr/J) 



(28) 



where C is some constant independent of magnetic field. 

If we put T = in Eq. l|^ we get the existing formula 
for the two-pulse echo amplitude Q. To illustrate this 
formula we plot on Fig.[4]the echo amplitude as a function 
of the external magnetic field for several values of the 
ratio p = T/r. For T = the magnetic field dependence 
is the same as for two-pulse echo (Fig. |4^). But with 
increasing of the ratio T/t the echo amplitude becomes 
larger (about factor of two) and gets more additional 
peaks as a function of magnetic field. 

As follows from Eq. (|28p the absolute maxima of the 
echo amplitude take place at 



2J 



kn 
TTp' 



fc== 0,1,2,.. 



(29) 



From positions of these maxima one can calculate mag- 
netic dipole moment fi of tunneling nuclei. Maxima at 
fiHT/2J = kn which exist also in a two-pulse echo signal 
are very shallow and can not be used for this purpose. 



VI. DIPOLE-DIPOLE INTERACTION 

Let us now consider the case where fine splitting of 
TLS levels is determined by dipole-dipole interaction of 
its nuclear spins 



E 



(1.2) 



(30) 
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Here fii = fiiJi/Ji magnetic moment operator of ith nu- 
clear, 3i is its spin operator, H is the external magnetic 
field. Radius vector = r[j'^\ from zth nuclear to jth 
nuclear depends on the well number of TLS (1 or 2). The 
summation in Eq. (|30[) is taken over all tunneling nuclei. 

The spin operators Wj and Vj in this case arc symmct- 
ric and antisymmetric parts of spin hamiltonian Ej ' 
correspondingly 



Wj 



E 



(1) 



E 



(2) 



Vj 



E 



(1) 



E 



(2) 



(31) 



Note that the first term of spin Hamiltonian ([50]) . which 
is associated with the external magnetic field (the Zee- 
man part), does not include vectors Vij that change under 
transition of the TLS from one well to another and, there- 
fore, does not contribute to the antisymmetric part of the 
Hamiltonian Vj. The dipole-dipole interaction (the sec- 
ond term of Hamiltonian ((30)) ) enters both operators Wj 
and Vj . Finally we should perform the unitary transfor- 
mation similar to the previous case, Eq. 1251 

To simplify the problem, let us consider further the 
limit of high magnetic fields when Zeeman energy is much 
bigger than the energy of dipole-dipole interaction Ed 
and Ed{T + T)/h < 1 which is similar with (gT]). In addi- 
tion we restrict ourself the case of dipole-dipole interac- 
tion of identical nuclei with spins 1/2 as it takes place in 
glycerol. Then by analogy with [10] we get from general 
formula, Eq. ([2T|) . for the echo amplitude the following 
expression 



by dipole-dipole interaction of TLS nuclear spins Vj. We 
see that Eq. ([5^ is formally identical to Eq. (^5]) for 
the case of quadrupole interaction with J — 1/2. This 
coincidence is due to that we have considered a dipole- 
dipole interaction of two identical nuclei. It disappears 
in a more general case. 



VII. CONCLUSIONS 

We developed a theory of the three-pulse electric- 
dipole echo in glasses in a magnetic field. As usually two- 
level systems contribute to the echo signal at low enough 
temperatures. We have shown that magnetic field depen- 
dence of the echo amplitude is due to quadrupole electric 
moments of tunneling nuclei and/or dipole-dipole inter- 
action of their nuclear spins. In the case when Zeeman 
energy is bigger than the quadrupole energy or energy of 
dipole-dipole interaction the magnetic field dependence 
of the echo amplitude appears to be a simple universal 
function of magnetic field independently of microscopic 
structure of TLS's. In comparison with the two-pulse 
echo the magnetic field dependence has a more reach 
spectrum of echo oscillations. From positions of the max- 
ima of echo amplitude as a function of magnetic field one 
can find the magnetic dipole moment of tunneling nuclei. 
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sin^ ifiHr) si-a^ {^iHiT + r)) 
sm^{2nHT) sm'^{2fiH{T + t)) 



(32) 
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The coefficient C in this expression is now determined 
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